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ABSTRACT 


A model is developed for relating continuous review 
inventory policies for repair parts to system availa- 
bility. The system consists of S identical unit systems, 
each of which is a series of k-out-of-n structures. Unit 
system states are zero or one. An optimal cannibalization 
policy is assumed. Under this assumption the number of unit 
systems up is always the maximum possible for any given 
vector of backorders for the N part types in the system. 
The distribution of backorders under a (Q,r) policy with 
Poisson demands for each part type is used to derive 
expressions for system availability as functions of the 
QO, £ vectors. For sampliciey at 1S) assumed. cae order 
quantities are set by an operating level in terms of days 
of supply. A numerical technique is presented for finding 
the vector of reorder points (safety levels) which mini- 
mizes the expected cost of on hand inventory subject to one 


of two alternative availability constraints. 
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Il. INTRODUCTION 


A. BACKGROUND 

In the management of logistic support for weapons 
systems a significant problem in resource allocation is 
the determination of optimal inventory policies for repair 
parts. Inadequate stocks of repair parts result in low 
system availability. Stocks sufficient to insure with 
high probability that any part will be immediately avail- 
able when needed my tie-up more resources than are 
justified. 

An indication of the concern of the logistics manager 
with this problem is the following statement from the 
revised General Objective No. 1 for the Navy Supply Systems 
Command: "For technical material, optimum support is that 
which minimizes downtime of weapons systems due to lack of 
repair parts and components." [Ref. 1] In order to accom- 
plish this objective, techniques are needed for relating 
repair parts inventory policies to weapons systems 
eave lability. 

The mathematical models designed to deal with the impact 
of repair parts shortages have usually been formulated as 
a minimization of total cost, where the total cost includes 
ordering, holding, and backorder costs. In theory the 
backorder costs are a measure of the impact of shortages on 
system availability. A difficulty in practice is that 


backorder costs are hard to estimate. One technique for 


ime 


finding a value of the shortage cost is to impute it based 


on a required maximum probability Pa of being out of 


ee 
stock. See Deemer and Hoekstra [{Ref. 2, p. 5] for a dis- 
cussion of this technique. The required value of Eid is 
determined by management judgement. This judgement pre- 
sumably includes some intuitive consideration of the effect 
of repair parts shortages on system availability. 

The total cost minimization formulation does not, how- 
ever, permit any direct correlation of inventory policies 
and system availability. Inventory policies computed using 
that formulation are determined for each part independent 
of all other parts. System availability is, however, a 
function of (Q,r) = (Qy eX y)reees (Qu Ty)- Therefore, if 
we wish to correlate inventory policies to system avail- 
ability, we need a multi-item model rather than one which 
deals with each item one at a time. 

There do exist a number of techniques for explicitly 
correlating system availability and single period inventory 
policies. One of these is the optimal redundancy approach 
discussed in Chapter 6 of Barlow and Prochan [Ref. 3]. The 
problem discussed there is one where at the beginning of 
the period there is a quantity a of each part type j 
on hand, and resupply is not available until the end of 
the period. The object is to determine the value of each 
n. so that system availability is maximized subject to 


some cost constraint or, alternatively, to minimize the 


cost of achieving a required availability. The cost of the 


eZ 


parts inventory may be interpreted as dollars, weight, 
volume or some other measure of the amount of a "resource" 


which is used up by the inventory. 


B. GENERAL NATURE OF THE MODEL 

The approach taken here is to apply the redundancy 
optimization ideas to a continuous review inventory sit- 
uation. In particular the model addresses the problem 
of optimal parts inventory policies in support of a weapons 
system at the direct support echelon. The objective 


function to be minimized is the expected cost of on hand 
ST TT ae sees ettitanea fen 


ee Sgn ane my Se ti 
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The model is an idealization of the repair parts 
supply support furnished by the direct support maintenance 
unit for the population of a particular major item in an 
operational Army unit. The system consists of S identical 
unit systems, which might be tanks, aircraft, howitzers, 


Or some other set of equipments of the same make and 


model. 


igs 


C. ASSUMPTIONS OF THE MODEL AND THEIR IMPLICATIONS 
Before turning to the formal description of the model 
we shall discuss some of the important assumptions and 
their implications here. 
i. The Inventory System is Ssinglestcnowen 
To simplify the problem we assume that the supported 
military unit does not carry a stock of repair parts, but 
that parts are immediately available from the direct support 
maintenance unit if there is stock on hand at the time a 
demand occurs. The inventory system thus is single echelon. 
2. Procurement Lead Time is Constant 
The procurement lead time (order and ship time) is 
assumed to be constant for each part type, where part type 
indicates a particular Federal Stock Number. Lead times may, 
however, be different for different part types. We further 
assume that the source of supply is never out of stock. 
Thus, the model ignores variability of lead times and the 
possibility of reducing lead times by using a higher requi- 
Si tl One r homey 
3. Unit SyvStemiscace: Valucs aren sc Lomeli. ne 
The supported system consists of S identical unit 
systems. Each unit system is assumed to exist in one of 
exactly two states depending upon the states of the in- 
stalled parts: either it is up, i.e., it is fully capable 
of operating satisfactorily; or it is down, i.e., totally 
ineffective. This assumption means that the possibility 


that the unit system may be partially effective is not 
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considered. Actual equipment does exist in partially 
effective states. For example, if the 50-caliber machine 
gun on a tank is not working, the capability of the tank 
is reduced in some of its roles, but it is obviously not 
totally ineffective if all other components of the tank 
are working. The zero-one assumption is made for mathemati- 
cal simplicity. The two state categorization is. not too 
different from the equipment serviceability code (ESC) 
ratings currently used within the Army for determining the 
materiel readiness of Army units. Possible ESC ratings are 
green, meaning fully operational and capable of operating 
in combat for 60 days; amber, meaning fully operational and 
capable of operating in combat for 30 days; and red, mean- 
ing not fully operational. We have simply dropped the amber 
category. 
4. Component Part State Values are Zero or One 

Similar to the zero-one assumption for each unit 
system two possible states are assumed for each installed 
part: either it is working or it has failed. Thus the 
part is 100 per cent effective or totally ineffective. 

5. The Unit System is a Series of k-out-of-n Structures 

As mentioned in Section 3 above, the state of a 
unit system is determined by the states ae ats component 
parts. If there are a, parts of type j initially 
installed on each unit system, it is assumed that there is 
a number b. equal to or less than a5 such that at least 


my, parts of type j must be working if the unit system is 
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to be up. If the number of working type j parts is less 
than Dia the unit system will be down. Further we assume 
that, if the number of working type j parts is equal to 
or greater than b. for all Jj, then the unit system will 
be up. These assumptions may be summarized by saying that 
the unit system 1s a series of k-out-of-n structures, as 
defined by Birnbaum, Esary and Saunders [Ref. 4, p. 58]. 
Oug D5 corresponds to these authors' k and our a. £6 
tee lyre. 


6. The Number of Unit Systems Up is the Minimum of the 
Number Upt@wilthm@respcct comedem mart Type 


If 2Z is the number of unit systems up at some 
arbitrary time and a is the number which would be up if 
all non-working parts other than type j were replaced 
by working parts, then Z is called the state of the 
system, and Zs is called the state of the system with 
respect to part type Jj. We assume that at any time t the 
state of the system Z equals the minimum over j of the 


fi cg 
J 


7. An Optimal Cannibalization Policy is Observed 
Perhaps the most significant assumption of the model 
developed here is that a policy of optimal cannibalization 
is observed. By optimal we mean that given any vector 
W = (Wo,.-- Wy), where = denotes the number of working 
parts of type j, after the cannibalization operation the 
maximum possible number of unit systems will be up. The 


mathematical structure of the supported system is an 


adaptation of the structure of the cannibalization model of 
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Hirsch, Meisner and Boll [Ref. 5]. The reason for using 
che S Strverere Fs that, "iver Gy weeter™ ¥ = (Yireee Yu), 
where Yi denotes the number of backorders of part type 

J, there is a unique number systems up. This would not be 
the case if cannibalization were not allowed. Consider, 
for example, a system in which each unit system has two 
parts of type j installed, and both of these must be 
working if the unit system is to be up. Suppose there are 
two backorders for this part type. This means that there 
are two non-working parts of type j in the system. If 
one of these were on one unit system and the other were 

on a different unit system, then the number of unit systems 
down for this part type would be two. If, however, both 
were on one unit system, then only one unit system would 
be down for this part type. Under a policy of optimal 
cannibalization there would always be only one unit system 
down for this backorder situation. Since the probabilities 
for the number of systems up used in this model are com- 
puted based on the probability distributions of the = 
the cannibalization assumption is an essential feature of 
the model. The degree to which cannibalization is actually 
used to increase system availability in practice depends 
upon command policy and the practicality of taking parts 
from one unit system to make another operational. If the 
supported unit is spread-out geographically, it may not be 
feasible to cannibalize. Also the maintenance effort in 


removing and replacing parts under the cannibalization 


Pal, 


operation may overload the capability of the maintenance 
Organization. Nevertheless, some cannibalization may be 
the best strategy for optimizing system availability when 
there are constraints on the quantities of repair parts 
the maintenance unit can carry. 

8. The Demand Distribution foresEach Part Type lis Pome 

The demand distribution for each part type is assumed 

to be Poisson with mean equal to the product of the failure 
rate for the part type times the number of parts of that 
type in the system. The failure rate of each part type is 
assumed to be the same in all of its applications. The 
Poisson assumption implies that the mean time between 
failures of a given part type is exponentially distributed. 
It is true that the exponential distribution may be a poor 
fit for the mean time between failureswhere an individual 
part is installed and replaced with a new part of the same 
type immediately upon failure, particularly when the part 
is subject to wearout. Cox [Ref. 6, p. 77] indicates heme 
ever, that the pooled output of a number of renewal processes 
tends to have the properties of a Poisson process as the 
number of renewal processes being pooled gets large. Con- 
sequently, the Poisson assumption may be a good approximation 
when the total number of parts of a given type in the 
system is large. Another implication of the assumption is 
that the distribution of demands is not affected by the 
number of unit systems down. This in turn implies that an 


inoperable unit system continues to generate part failures 
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even after it has gone down. Errors introduced by this 
assumption will not be too great if the availability rate 
is reasonably high or if failure rates are estimated on 
demand histories over a period when the availability rate 
was about the same as the required availability rate for 
the period for which inventory policies are being computed. 

9» Demand Rates are Linear with Activity Levels 

The parameter in the demand distribution is assumed 

to be a time rate. As discussed by Soland [Ref. 7, p. 45], 
the natural parameter of the demand process may be miles 
driven, rounds fired, or some other measure of the usage 
or activity level of the system. The device used in the 
model to account for this fact is an "activity level" 
multiplier L. for the pooled demand rates for each part 
type. Suppose, for example that rounds fired were the nat- 
ural parameter for the demand presse Further suppose 
that a certain number of rounds were fired on the average 
by each unit system per month during the period when 
demand data were accumulated for estimating failure rates. 
Now suppose that during the period for which inventory 
policies are being computed that the programmed number of 
rounds per unit system per month is doubled. Then the La 
would be two, and the failure rate estimates for each part 
type would be doubled. Use of this device assumes that 
failure rates are linear with the level of usage. This 
assumption is not in general true, but any more realistic 


means of handling this problem would complicate the model. 
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For part types not subject to aging, this assumption is 
probably fairly good. 
10. Parts are Replaced Immedtavelyeaupon Pal lure 

This model assumes that each part is replaced 
immediately upon failure with zero replacement time if 
a spare part is on hand at the time of failure. Also, if 
backorders exist for a particular part type, the failed 
part will be installed immediately upon receipt by the 
maintenance unit. Of course, these assumptions are not 
true in the "real world." To account for the fact that 
repair times are finite the following procedure could be 
used. First estimate the availability which we would 
expect to achieve if there never were any shortage of repair 
parts. Call this estimate the availability with respect 
to maintenance. Call the estimated availability predicted 
by this model the availability with respect to supply 
Support. Then by taking the product of these two estimates 
we Should have a fair estimate of the overall availability 
to be expected given the inventory policies and the repair 
capabilities of the maintenance unit. In any case, the 
reader should keep in mind that availability values computed 
with this model need to be corrected for down time due to 


non-zero repair times. 


20 


It. MATHEMATICAL MODEL 


~. FORMAL DESCRIPTION OF THE SYSO Rem STRUCTURE 

The system consists of S identical unit systems. On 
each unit system k there are n loci where parts are 
installed. Let the ane denote the ith locus on the kth 


unit system and let the set 


Nn 


} 


i=1 


A, = fh 


k Kea 


denote the set of loci on the kth unit system. Let the set 


ie 
A=) oj A 


k= fs 


ft Ne woe 


denote the set of all loci on the S unit systems. Let 
Xs be a random variable indicating the state of the ith 
locus on the kth unit system, such that 


Xe 1 Liflocus AR ;, ~copmbalmse~aowerking part of 
the type required in that locus, 


0 if locus Ayy fails to contain a working 
part of the required type. 


Now = (Xpress XL) 1s a vector of zeros and ones which 
describes the state of the kth unit system. The possible 
values of Xy correspond to the vertices of the unit cube 


in Euclidan n-space. Let x1 denote the set of all pos- 


nN 


eile values of eel ze GL is 2. Let the 


k Ak 


random vector 
aS (X)7-++ 1X) 


ae 


be an nS-component vector indicating the state of the 


entire system. Let », denote the set of all possible 


values of 2@apeThems2 zouot® as Des Let the set 
n 
¢ = {FSN.} 
z) j=1 


represent the set of different part types installed in the 

n loci of each unit system. Further let G(FSN.) be the 
subset of AL. which has as elements all those loci requir- 
ing part type pon It is assumed that if locus Ay 4 ©GSN}) 


f 


then ee Sy 1 aif and only if a working part of type —? 
is installed in that locus. This assumption implies that 
different FSN's are not substitutes for one another. Let 
the size of GtEeNTy = au 1.e., the number of applications 
of FON. on each unit system is a Assume that for each 


j there exists a number De such that 


ile ese Danese a: 
J J 

2. If the number of working parts of type et in- 
stalled on unit system k is equal to or greater 
than Das the unit system will not fail due to 
PON 

a] 

3. If the number of working parts of type FSN. Loi 

stalled on unit system k is less than se the 


unit system will fail due to ee and 


4, oe 1s independent of bs for@ali. 27 7: 


teu, che State Of the Unie system with Mes UeCiE LO FSN4, 


denoted OK? be defined as the state in which all loci 
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net Gneehe set Seer have working parts installed. 

In other words ies describes the state of the unit 
system if ne were only part type subject to failure. 
Tf ms is a random variable indicating the number of 


working parts of type — installed on unit system k, 


then 


U5 = lif Wy Zs Das 


0 otherwise. 


Nem@ let the state of the isvstemewith mespect te Se 
5 


denoted ary be defined as Lae Z. indicates the num 
ber of unit systems which would be operational for a given 
set of values of eae k = l,.--,S i1f the only fallable 
part type were See 

For a given value of the system state vector xX it 
may be possible to increase the value of 2 by cannibali- 
Zation.. We assume that for any value of X cannibalization 
will be effected in such a way that the maximum value of Z 
mer that X will be achieved. See Hirsch, Miesner and Boll 
[Ref. 5, p. 336-342] for a detailed description of the 
cannibalization operation for the type of structure being 
discussed here. Let 2Z* denote the state of the system 
after an optimal cannibalization, and let as denote the 
state of the system.with respect to i, under a policy 
of optimal cannibalization. 


Note that oe may be considered to be a function of 


tae the number of working parts of type j, since a is 


aS 


the value of the number of operational unit systems after 
an optimal cannibalization when the only part failures are 
at those loci associated with part type j. It is assumed 
that Zeki) is a monotone non-decreasing function of Ws: 
Let at 2} denote the set of values of the number of 
working parts ne such that the value of = is at least 


as great aS Z, 1.€e., 
K(W.,z) = {W.: 0 < W. < Sa. and Z2*(W.) > 2}. 
( ; z) { ; Sellen a. an 4 ;) z} 


Note that 


K(W, OVC K(W.,1)C. : -CK(W,,S)C K (W, ,S+1) , j=l aa 


and since max a S, 


iI 
= 


me 


The largest value of z such that a ois S. 
Further note that the minimum value of Ws such that Z 
ae =z is 94 since 2s is the minimum number of part 
type j needed on each unit system if that unit system is 


to be operational, in symbols 
7b. = mim {wees Woe kere) eee 
J J J J 


Let Kex be an arbitrary set where x is the set of all 


system state vectors X. The simetea tor Stoner ren I, 1s 


defaned for all K and X as follows 


I,(X) = 1, if Xex, 


0; “TE xen & 
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Now let us use the symbol tw. 2 zb,] to denote the set of 
aan? eX "Steh “that Ww. 3 zb., i.e., (Keg: w. (X) = zb.}. 


Hirsch, Meisner and Boll show (using somewhat different 


notation) that if 2Z2* = min Z*, then 
J 
S S N 
z~ = f z I I ; 
k= fw. > zb ,} 


N 7 
1 7 fw.2zb.} k=l 4j=1 
7 


and the authors call this the "representation theorem" [Ref. 


by ©. S89]. If S 


1, 1.e., we have only one unit system 
then Z2* = U, and U= 1 if and only if each one of the 


indicators Ty is equal to 1. The structure of 


w.=b.} 
the unit system in this model is thus a coherent structure 
in the definition of that term given by Birnbaum, Esary and 
Saunders [Ref. 4, p. 61]. Further each unit system may 

be considered to be a k-out-of-n structure with respect 

to each part type j, where a k-out-of-n structure is one 
which has n parts of a single type and is operational if 
and only if at least k out of the n parts are working. 
Also, we can consider each unit system to be a series system 


composed of N_ k-out-of-n structures, where the k's are 


the Paes and the n‘'s~ are the py Bx 


B. CORRELATION OF INVENTORY POLICIES AND SYSTEM AVAILABILITY 
Consider now the effect of inventory policies on the 

expected value of the number of operational unit systems 

E(Z). If there were an infinite number of repair parts of 


each type available, the value of X would always be 
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1 = (1,1,...,¥), since it is assumed thateach parce 
replaced immediately upon failure as long as spares are 
available, and the state of the system 2Z would always 

be S. There are in fact various constraints on the ability 
of the support unit to carry inventories of parts. The 
result is that from time to time demands will occur for 
which no replacement part is immediately available. For 

any state vector xX, dene@te the stares Scere gso. roma 
optimal cannibalization as xX* and the value of 7% diem 


Cannibali zakion as ~2* ~ 7 Then 


XAE ) = (XT 7 ft). XB CE), ic ey U 


—_— g 


is a stochastic process where for each fixed t, Xs is 
a random variable taking the values zero or one. Now for 
each subset «K s@f mm, the set of pessabie valmesmot ae 

let PB, (K) denote the probability that 2@(t) iIsianpaie 


ment of the set K, lLe., 


Pi (k) = Prix*(t)eK}, KCx . 


As shown by Hirsch, Meisner and Boll the probability 
distribution ‘defined above concentrates all of its mass on 


the set of maximum points M which are the possible 


we 
vectors X*, the state vector after an optimal cannibali- 
zation. Any vector xX in which the number of working 
parts of each type is given by the vector W can be 


transformed into another vector x! which has the same 


value of W. Thus, since an optimal cannibalization yields 
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_——_ —————————— 


Le 


the maximum value of the system structure function Z for 
eeciiven cartce VEStTOM 2G, Sthere is a timmqre valuerot the 
state of the system after an optimal cannibalization 2Z* 
fee each valuewof the wector W = (Wir---/Wy). The set of 
all X corresponding to avgiven W vector is sealled an 
equivalence class, and the set My. 1s a subset of this 
equivalence class where Z asa Susteren Of "* “takes” its 


maximum values over the equivalence class. We denote each 


~& in My as, X*. Te d@scravestite variation in time of 


the system state function 2Z, set 
Zea) = ones (e)) cee OF 


and we note that the probability distribution of 2Z*(t) is 


given by 


Pr{zZ*(t) = z} Pr{Z*(x*(t)) 2 az} 


P,{X: Z*(X) 2 2} . 


Let ee) represent the number of working spares of 


type j at time t. Under the assumption that Z* = min 2*, 
J J 


we have immediately from the representation theorem that 


S N 


d 


S 
Zaye LT = ee : 
- {W,(t)2kb.} ‘een (W(t) 2kb 5} 


N 
ia 
( 
fa 


Assume that 


i. Bor each index .j = 1,...,N, the loel tewhich 


parts of type pee occur are indistinguishable in their 


Zl 


effects on the lifetimes of the parts installed in them, 
i.e., the failure rate at each instant of a given part type 
ae does not depend on the locus in Sauer in whaleh 
the part is installed, nor on the particular sequence of 
loci through which it has passed. 

2. Parts operate independently, i.e., the lifetime of 
a given part is not related to the lifetimes of any other 
parts. These assumptions make it reasonable to postulate 
that the joint distribution of (W, (t),..- Wy (t) ) does not 
depend upon the particular cannibalizations involved in the 
process {Z*(t), t > 0 and that the N_ stochastic pro- 


ces cto 
{W, (t), Leon, Pei {Wy (t), aon, 


are mutually independent. Thus the expected number of unit 
systems up under optimal cannibalization at time t, 


E(Z*(t)), is given by 


Ss N 
Bia * (C)) ee een ] 
1 {W.(t) 2 kb.} 
S N 
= © I Pri{W.(t) 2 kb.}. 
i= 9a J : 


Further, the probability that Z* (@)equatls jaegleaa =z 


1S given by 


Priz*(t) =z} = 


Ia 


Priw. (t) oo zb.}. 


j=1 
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Now let co be a random variable denoting the number 
of backorders for part type Pom at time t. Under the 
assumption that every failed part is replaced immediately 


by a working part from the inventory of spares, we see that 
ree) — See Le 
5‘ ) 5 5‘ ) 
and 


{W. (t) Ps kb, } = {Sa,- Y(t) 2 kb.} = {¥. (t) & Sa,- kb}. 


Hence 
N 
Prizg*(t) > kK} = a Prty tt) 2 Sanne . } 
4=1 J J J 
and 
S N 
BE(Z*(t)) = £ eo Pay . tt) >See —» kis } 
Ke jad J J J 


We assume that the inventory system is continuous review 
with a (O.,r4) policy for each FSN, where Q. denotes 
order quantity and ry denotes reorder point. We also 
assume the lead time demand is Poisson distributed. Hadley 
ema Whitin [Ref. 8, p. 184} show that un@er a (Q7r) “policy 
when the lead time demand is Poisson with parameter oT 
that the probability that there are y backorders at an 
arbitrary time t under steady conditions may be derived 
me follows. The inventory position, IP, is défined asthe 
stock on hand plus stock on order minus backorders. The 
IP varies between r+l and r+Q, where r is defined in 


terms of the IP. Further, Hadley and Whitin show that the 
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probability that ‘the IP TYsS"i'n” any™staee wrt], ate eee 
is 1/Q. Thus if the inventory system is in state r+j at 
time t-t, the probability of y backorders at time t 
equals the probability that ytr+j demands occur during 
the period (t-T) to t for y>0. The probability that 
y=0 is the probability that the demand during (t-T,t) 

is less than or equal to r+j. Hadley and Whitin do not 
make this distinction for y=0; consequently, the formula 
they derive is not valid for y=0. When 0O is substituted 
for y in the expression they derive the resulting value 
is the probability that y=0 and on hand inventory = 0. 
Terms for the probability that y=0 and on hand inventory 
is greater than zero are left out of their expression. The 


expression Hadley and Whitin derive is 


QO Vt +o 
p= (1/0) £ plytrt+jret) = (1/2) 2& pila om) 
j=l u=ytrt+l 
= (170) [P(yv¥relyot) — P(yt+rt+Gri;et) 1, 
where p(u;opTtT) is the Poisson probability mass function with 
parameter. ot at the point u for uu = 0,1,2,..2: aga 
P(u;ot) = 2% plvret) is the Poisson complementary cummulative 
V=u 


distribution function. The above expression is valid y 


Caeateiawthan 27650. FOR vo — 09 
QO Sr) 
mii) = (ZO. ta EG ion) 
j=1 u=0 
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give) rene tet tale te) | 


il m0 
- 


J 


Q 
eo) ee 


P(r+j+1:pT). 
1 


a 


We assume that the pooled demand rate 5 for all 


parts of type j is given by 


where Lo is the activity level multiplier discussed in 
Section I, S is the number of unit systems, a. is the 
number of applications of part type j on each unit system, 
and p. is the failure rate for part type j when L 


a 
"Ss one. 


Let Ts be the procurement lead time for part type j. 


Demand is assumed to be Poisson distributed with parameter 


ftr.. Thds in term if ; ; 
p=. S in temere va lya), 

sa. ~kb 5 

PRiz<(t) > P= 2 Woy.) . 
J y .=0 7 2 
J 
Also since 2*(t) = min aaa and the i are assumed 
J 


to be independent, the probability that Z*(t) is equal 
to or greater than k is the probability that each a 


is equal to or greater than k, i.e., 


a 
Prw*(t) > kK} =" TI y 
J 


a2 


Now the expected value of a non-negative integer random 
variable is the sum from 1 to ~» of the complementary 
cumulative probability function; thus 

S 


{Pr Z*(t)e > khe= DeboaeZ* (tc) See 
i eal 


Bae) se 
k 


lil ws 8 


Since the probability that 2Z2*(t) is greater than S is 
zero. Therefore, the expected value of the number of 
operational systems in terms of the probability mass 
functions of the backorders for each part type is: 

S N Sa.-kb. 


E[z*(t)] = ©. nm 2 2? 


Walys), 
kel j-1 y.=0 Jd 


where: 


a 
J. 
ke 
J. 
ll 


P +x .t - Pp +tr.t+0O., : 
Cro) 1 ty os Ll) (y. t. O,+1)], for yi70, 


ana 


Q 
. 12> @ £ P(r.tmtl), £ _g 0. 
eT (1/Q) a (eye); foray. 0 


— 
he 
II 


The een of optimum inventory policies may be formulated 
as the minimization of the expected cost of on hand inven- 
tory subject to one of the following constraints: first, 
that the expected number of operating systems E(Z) must 
be equal to or greater than some required fraction of the 
total number,or, second, that the probability that at least 


k of the unit systems are operational must be equal to or 
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greater than a minimum assurance level. The expressions 
used to relate inventory policies and these constraints 
are those developed above for E(Z) and Pr{Z2k}. Note 
that since the expressions refer to "steady state" con- 
ditions, and since optimal cannibalization is always 
performed, the symbol 2Z will be used henceforth in lieu 
of Z2*(t) and Z will indicate the "steady state" random 
variable 2Z*(t). 

It turns out that for part types with low demand rates 
the optimal policy under the formulation of the problem 
indicated above is not to stock these low demand items 
at the direct support level, but to order them from the 
source of supply as demands occur. If this is the case, 

r will be set to -l and Q tol. This means that an 
order is placed as soon as a failure occurs. [In this 
Situation the distribution of “ reduces to the Poisson 
distribution with parameter PaTss Since the inventory 
position is always equal to zero, and there is, therefore, 
only one state for the inventory position. Hence the 
probability that there are exactly Ys backorders at any 
arbitrary time is Ea rec ee 

It is interesting to consider what would happen if 


0 were small for all N part types, and no stock were 


if. 
7 os 
carried at the direct support level. Assume that each part 
type were always available from the source of supply with 


a procurement lead time of Tay a constant, which might be 


different for different part types. Now, if no parts are 
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carried at direct support level, but are ordered only when 


demands occur, 1.e., t. = -l, Q. = 1 for all j, then 
S Sa j-kb, 
Pr{Z >k} = I bi Dil Ves OniTaiy) 
j=1 y.=0 anes) 
J 
and 


S N Sa4-kb., 
ey: J 
k=l j=l y5=0 


ll 


E(Z Oe lee) es 

(Z) ony p a 

The values of the above functions can be used to estimate 
the lower bounds for E(Z) and Pr{Z >k}, when demand 


rates and procurement lead times are Py ey j=l. 0 3 pe 
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III. OPTIMIZATION TECHNIQUE AND NUMERICAL ANALYSIS 


-m FTE COST PUNCTION 
The objective function to be minimized is a cost 


function of the following form: 


N 
Bic) = 2F C.BWO.7r.), 
j=1 x — 
where 
E (C) = Expected cost of on hand inventory in 
dollars, cubic feet, pounds or some 
other measure of a resource in short 
supply, 
os = The unit cost of part type jj, 
D(Q. rs) = The expected on hand inventory of part 


type j, given order quantity Q 


and reorder point Ea 


J 


If E(C) is in dollars, the above function gives the 
expected amount of funds tied up in repair parts inventory 
for the system at the direct support echelon. If E(C) 
were in pounds the function might represent the expected 
amount of load carrying capacity of the maintenance used up 
by the repair parts for the system in question. 


The value of DME) 1s given by 


-7,0.) = . = Pay B(Q.,r.), 
D(Q. 4) sons lial 4 vr. Pats + (Q. v3) 
where 
B a. = -£ -Y. -7Sds 70s 
(Q. = 4 Y «Waly Q. r5) 
J 


be 


is the expected number of backorders, given the inventory 
policy (Q.,r5). Hadley and Whitin [Ref. 8, p. 184] show 


that 


B(Q,r) = (1/Q) = y[P(y+r+l;ot) - P(y+rt+Q+l1;pT)]). 
y= 0 
Bus OPTIMIZATION TECHNIQUE 


1. Minimizing Cost Subject to a Required Expected 
Number of Unit Systems Up 


Suppose it is required that the expected number of 
unit systems up must be equal to or greater than AS, 
where A. is the required availability with respect to 
supply support, 0 < A. <1, and S is the number of unit 
systems. Further, suppose it is desired to minimize the 
cost of expected on hand inventory needed to achieve 


E(Z) 2 AS. The problem may then be stated as 


N 
Minimize ew | Coli. | Cue) 
(L770... Ln) j=l J J J 
(Qy eee 1 Qy) 


subject to # ©b(Z) 2 ASS. 


Finding the optimal solution to this problem for 


(Qarta)s j =1,...,N is a formidable task due to the com- 


plexity of the cost function and “the [funetven for the 


expected value of 2Z in terms of the Q. and es: Te 


Simplify the problem we have chosen to set the values of 
Q. by an operating level of supply in terms of days or 
months of supply. Thus we set 
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Q. = UP ole, +" Ey pes tosh, ~ D5 >A 


= 1, otherwise, 


where: [(u] = The greatest integer equal to or less than u, 
L, = The operating level of supply, e.g., one month 
of supply. 
The problem now becomes: given the values of Q. determined 
as indicated above, find the vector (Lp pe++ Ky) which 


Minimizes the cost function subject to the required level 
of availability. We shall restrict the er to be: -»- FT: 
After computing the Q. values, the next step in the pro- 
cedure is to compute a table of values of V5 lyai r5 = -l) 
for all Jj. This table is useful in that the values of 
ty; aw may be found for other values of = as 
follows: 


ratl 


s QO; ‘ = ), ‘ ? -l , £ 2 = OQ, 
W007 Fs) J oe a2 te Y. 


os 

< 
rs 
i 


+tr.+1; -l), for w> Ow 
WAY x. ) Y 5 


The number of values needed for each V5 depend upon the 
demand rate and the accuracy desired in the computations. 
After the table of yw values is computed the next step 1s 
to find a set (rine. .oud) of t. such that Eiv@.,}, EME 
expected value of 2Z with respect to part type j, is 
equal to or greater than A_S, and such that Pr(Z.=S;r,) 


J 
1s concave for ts = rie Each eS must be large enough 
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to insure a = A_S since ENeny gives the expected 
number of operating systems when every other type except j 


is always in supply. E(Z) can not be greater than Bae) 


because it 1s assumed that at any arbitrary time =min a4. 
J 


Now ea. is given by 


S Saj-kb, 
| Ac vk Vey. ee e. 
J ve i el y 
J 


Using this expression find the minimum - such that 


ae 2 AS. (See Appendix A for a discussion of the need 


aa J 
After the set of initial values of the r is found, a 


to require each Pr(Z. > k; oe to be concave for Ee 


marginal analysis technique similar to the approach of Karr 
and Geisler [Ref. 9] and to one form of optimal redundancy 
algorithm of Barlow and Prochan [Ref. 3, p. 166] is used to 
find at each iteration the part type which yields the 
greatest increase in E(Z) for the increase in the cost of 
the expected on hand inventory for that part type. On the 


first iteration! the Matic 


E(Z:r.t+l) - E(Z3r.) 
| © GMD alr fi) =D Gon 
; 5 bd, tr 5th) 3 {5 


is computed for each part type, where 


ACE ee ie 


Pr (22k) Pr(Z.2k; a 


us 


lun 


E(Z;r.+l) = 
J k 
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The maximum over jj of ie is determined and the ey 
corresponding to that maximum is increased by one. E(2Z) 
is set equal to Balhae ad coBpmespondangGuitO® th iis. Jeeeadi 
the new value of E(Z) is equal to or greater than AS, 
the procedure stops; if not, a new set of 4. is computed 
and the x corresponding to the maximum over j of the 
new a is increased by one. The procedure continues in 
this fashion until the value of E(Z) is equal to or 


greater than AS. 


2. Minimizing Cost SUbgect Oa Becgiumos Probably, 
that at Least k Unit Systems are Up 


Suppese that it iismmequimed. cham asP 1m eek 
P nin Ke and that it is desired to minimize the cost of 
the expected on hand inventory needed to achieve this 
probability. The general procedure for finding the optimum 
vector (Ly eee Ky) is the same as the E(Z) 2 ALS form 


of the constraint with the following expections: 


a. The initial = values (r¥,...,/rH) are the 
minimum values such that Pr(Z. 2k) > P pe (1) and 
Die ite alu Keae) gltSeeGONCAVCc DDT » qutO Gee > ot * 
J _—) J J J 


b. At each iteration be is computed as 


Pr{Z > K;Y.+l}—-DPYria 2 i Ge 
a C.(D(r.+1) - D(r. : 
J 5 = ) (rey) 


As in the first form of the constraint the r associated 
with the maximum As is increased by one at each iteration 


until DeahZe> kK)i> P 


> Pi, (Kiem Sance Pr{zZ 2k} is just one 


3 


term of the sum over k in the expression for E(Z) the 
computational effort at each iteration is much reduced for 
this form of the constraint, particularaly if the number 


of unit systems S is large. 


C. NUMERICAL ANALYSIS AND DEVELOPMENT OF FORTRAN PROGRAMS 
Two FORTRAN programs, ane for each form of the constraint 
have been developed and run for several sample problems on 
the Naval Postgraduate School's IBM 360 computer. The pro- 
grams are very Sr as and both use identical versions of 
the following subprograms: PSITAB, which computes a table 
of values of an area =Seaseetor j= 1,...5Ne7> PRZJK, a 
function subprogram which computes eee > Se EBO, a 
Function which computes the expected number of backorders 
for part type j give that the reorder point is tai 
PPT, a function which computes individual terms of the 
Poisson orSbabiaeey mass function; and PCUMT, which computes 
complementary cumulative terms of the Poisson distribution. 
With the IBM 360 it was found necessary to do the 
arithmetic in double precision because the single precision 
round-off errors for repeated multiplications were excessive. 
The table of Vig MY GE = -l) values was, however, stored 
aS a Single precision array to conserve storage space. When 
values from this table were needed they were converted to 
double precision with the standard function DBLE. 


Subroutine PSITAB stores the computed values of 


aS ee = -1l) in an array called PI(J,K), where 


40 





Pina, K) "= ¥z(k-lirs =). 


The array size is (N,100), where N is the number of 
different part types; however, for each J only values 
of PI(J,K) for K less than or equal to KMAX(J) are 


computed, where 


KMAX(J) = min J100, largest K such that 
PI(J,K) < 10714 


PI(J,K) for K greater than KMAX(J) are set to zero. 
For part types with low demand rates the number of values 
of PI(J,K) computed is much less than 100. The array 
Size used in the programs permits handling demand rates 
up to about 50 without significant truncation of the distri- 
Putien of oe 
1.» «Program EZMIN 

Program EZMIN computes the optimum Xr jwectoresmb- 
ject to a required minimum expected value of the number 
unit systems up. The main program reads the input data, 
writes the system parameters and the required availability, 
and calls subroutines PSITAB, INITAL, and OPTIMR. Sub- 
routine PSITAB computes the PI(J,K) array as discussed 


above. Subroutines INITAL computes an initial set of - 


values, rf,.--,ry, each of which satisfies the following 
conditions: co = Sir.) a concave function of & for 
i. > Ye; 
, es 
S 
QZ.) =—2 Pir@ieeewloneaye > Bes. 
y red J 
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Subroutine OPTIMR uses the marginal analysis technique 
previously discussed to find the Op Ua igeey cecum 
2. Program PKSMIN 
Program PKSMIN computes the optimum r_ vector 
subject to a required assurance that at least k out of 
S unit systems are operational. This program differs from 
program EZMIN only in the following respects. PKSMIN Main 
calls subroutine INTLZ instead of INITAL and subroutine 
OPTMZ instead of OPTIMR. Subroutine INTLZ finds an 
inteia set of rs values Cores TH such that Pr (2, 2k3r,) 
is concave for a“ 2 rs and such that = > a’ > P in’ 
Subrowtine OPTMZ computes the optimum Yr vector subjecuses 
the required assurance that 2Z be equal to or greater than 
ges 
3. Kkew Chartres 
Flow charts for program EZMIN are presented in 


Appendix B. Flow charts for program PKSMIN are presented 


in Appendix C. 
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IV. EXAMPLE PROBLEMS AND DISCUSSION OF RESULTS 


A. DESCRIPTION OF THE EXAMPLE SYSTEM 

The example system consists of 50 106-mm recoilless 
rifles, model M40Al. Several factors made this system a 
convenient example. First, the total number of repair 
parts is small enough so that the FORTRAN programs developed 
for solution of the two alternate formulations of the prob- 
lem can be used without modification on the IBM 360 G-level 
computer available at the Naval Postgraduate School. 

Second, Department of the Army Technical Manual 9-1015-221-35 
[Ref. 11] contains data from which failure rates may be 
estimated. Third, actual demand data for a six month 

period for 62 weapons used in training heavy weapons infan- 
trymen at Fort Ord, California, were made available to the 
author by the Fort Ord Post Maintenance Section. Finally, 
price data were available from the microfilm Army Master 
Data File Selected Management Data File [Ref. ll]. 

The parts list in TM 9-1050-221-35 contains 287 different 
Federal Stock Numbers (FSN). Of these only the 159 FSN 
which were indicated as. being combat essential were used 
for the sample problems. By limiting consideration to 
only combat essential parts we insure that some of the 
assumptions of the model are more nearly satisfied than if 
all part types were considered. For example, the model 
assumes that for each part type j there is a positive 


number Ds which is the minimum number of working parts of 
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type Jj needed for the unit system to be capable of per- 
forming its function. A part would not be combat essential 
fete m. were zero. Thus the assumption that the unit system 
is a series of k-out-of-n structures and that the possible 
states are zero or one is more nearly true if the unit 
system consists only of combat essential parts. 

The Federal Stock Number, the number of applications 
per unit system, and data for estimating failure rates 
were obtained from columns 2, 4, and 6, respectively of 
Section 2, Appendix B, Ref. 10. 
B. SAMPLE PROBLEM 1 - REQUIRED EXPECTED NUMBER OF UNIT 

SuolTEMSeUP GONSTRAINT 

The first sample problem was computed for the required 
expected number of unit systems up form of the constraint. 
Table I lists a summary of the input data and results for 


this problem. 


TABLE I 


SAMPLE PROBLEM 1 - SUMMARY OF INPUT DATA AND COMPUTED RESULTS 


NUMBER OF UNIT SYSTEMS 510 


NUMBERS OF SPARTA yPES 1S) 
TNE Ue REQUIRED AVAILABILITY 0.35 


DATA OPERATING LEVEL Vee 


ACTEVITY LEVEL 1.0 

EXPECTED NUMBER OF 

OPERATIONAL UNIT SYSTEMS 47.58 
COMPUTED EXEECTEDSCOSTAGE 

ON HAND INVENTORY $418.04 


RESULTS 
NUMBER OF ITERATIONS OF 
MARGINAL ANALYSIS PROCEDURE Pa 215 
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The above data show that for the system consisting of 
50 weapons with 159 part types and a required availability 
of 0.95 program EZMIN computed an r vector with an 
associated expected cost of on hand inventory of $418.04. 
The computed value of E(Z) is 47.58, which indicates that 
the marginal analysis technique overshot the required 
expected number of systems up by .08, since the required 
expected number was 50 x 0.95 = 47.50. 

Table II lists the input parameters and computed 
results for each FSN. Parts are listed in Federal Item 
Identification Number (FIIN) sequence, because they are 
listed this way on the Army Master Data File [Ref. 11]. 
The FIIN is the last seven digits of the Federal Stock 
Number. The data cards were sorted into this sequence to 
facilitate determination of unit prices. The columns of 
Table II contain the following: j, the sequence number; 
FSN, the Federal Stock Number; A= a., the number of 


J 


applications of part type j on each unit system; B= Dy, 


the minimum number of part type j needed for the unit 


system to be up; C=C., the unit cost of part type j; 


5 
RHOHAT = Par the failure rate for part type j; RHO = p = 
Sa.05, the pooled demand rate for all parts of type j in 


the system; TAU = Re the procurement lead time for part 
type j; Q= Qs, the order guantity; R= eof the reorder 
point; EB = 2 OE a the expected backorders for part type 
j; and EOH = D{Q. rs), the expected on hand inventory for 


part type a 
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C. SAMPLE PROBLEM 2 - REQUIRED ASSURANCE AT LEAST K UNIT 
SYSTEMS WE BCONSTRAINT 


The second sample problem was computed for the required 
assurance at least k unit systems are up form of the 
constraint. Table III lists a summary of the input data 


and results for this problem. 


TABLE III 


SAMPLE PROBLEM 2 - SUMMARY OF INPUT DATA AND COMPUTED RESULTS 


NUMBERWOF UNIT SYSTEMS 


50 
REOULRED NUMBER OF UNIT 47 
SYSTEMS UP 

REQUIRED ASSURANCE + opel 


OPERATING LEVEL 
ACTIVITY LEVEL 


PROBABILITY AT LEAST K O.29n 
UNIT SYSTEMS ARE UP 
EXPECTED NUMBER OF 47.46 
COMPUTED UNIT SYSTEMS UP 
RESULTS EXPECEED COST. OF ON $387.88 
HAND INVENTORY 


NUMBER OF ITERATIONS OF 
MARGINAL ANALYSIS PROCEDURE 





The above data show that for the same system considered 
in Problem 1 with a required assurance of 0.90 that at least 
47 unit systems are up program PKSMIN computed an r vector 
with an associated expected cost of on hand inventory of 
$387.88. The program overshot the required probability 
by 0.01 since the computed probability is 0.91 and the re- 


quirement was 0.90...After the. r vector which met the 


ae 


constraint was computed the program computed the resulting 
value of E(Z), the expected number of unit systems up. 

It 1s interesting to note that the expected number of unit 
systems up very nearly meets the requirement of Problem l, 
and the expected cost of on hand inventory is also about 
the same as for Problem 1. The computer execution time for 
Problem 2 was about 19 seconds compared to slightly over 
four minutes for Problem l. 

Table IV lists the input paramaters and computed results 
for each FSN. The column headings are the same as for 
Problem l. 

D. EXPECTED "COST OF ON HAND INVENTORY AS’ A FUNCTIGH Or 

REQUIRED AVAILABILITY 


1. Expected Cost Versus Required Expected Number of 
Unit Systems Up 


Program EZMIN was run for ten values of A. = Won, 
».-,0.99 in order to develop a functional relationship 
between the first form of the availability requirement and 
the expected cost of on hand inventory. Table V summaries 
the results. In addition to the expected cost of on hand 
inventory the table lists the number of iterations of the 
marginal analysis procedure required to find a solution and 
the number of reorder points which were less than the mean 
lead time demand for each level of A.- A reorder point 
less than the mean lead time demand is equivalent to a 
negative safety level. For many of these the computed order 


quantity is one and the computed reorder point is minus one. 
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This (Q,r) implies that the part is not stocked by the 
maintenance unit but is ordered from the source of supply 
as demands occur. 

The data show that expected costs of on hand inventory 
rise very steeply as the required level of A. approaches 


Uriel 


TABLE V 


EXPECTED COST, OF ONSEAND ZNYVENTORY VERSUS REGUTRED Ag 


Nexo f 
Negative 
A ; No. of Safety 
S EB (ec ) Iterations Levels 
O70 $16.80 oo 156 
0 oT > 22 602 64 Jab 
0132 $38.14 67 143 
Ome s Sue OT 140 T2 
0.94 $168.64 141 104 
Os $418.04 ZS ve 
05916 $645.76 238 Al 
Oro Sue22 a3 250 19 
0.98 $1619.00 322 ilu 
O32 $2406.83 386 2 


Figure 1 plots the expected cost of on hand inventory 
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Figure 1. Expected Cost of on Hand Inventory E(C) 
as a Function of Reguired Expected Fraction of 
Unit Systems Up 
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2. Expected Cost of On Hand Inventory Versus Required 
Number of Unit Systems Up with an Assurance of 90 
Pe meer t 
Program PKSMIN was run for six values of k = 45, 
»--,50 at a required assurance of 90 per cent in order to 
develop a functional relationship between the required value 


of k and the expected cost of on hand inventory at this 


assurance level. Table VI summaries the recsulice. 


TABLE VI 


EXPECTED COST OF ON HAND INVENTORY VERSUS REQUIRED 
NUMBER UNIT SYSTEMS UP WITH 90 PER CENT ASSURANCE 


Now*weF 
Negative 
No. of Safety 
Ik EGe) Iterations Levels 
45 Sen 108 150 
46 $84.97 102 i122 
47 $387.88 ale 85 
48 Seki, 4 7 218 49 
49 5 a OZ 278 let 
50 S3a19.92 S302 0 


Figure 2 plots the expected cost of on hand inventory 
E(C) as a function of the required number k of unit 


systems up with 90 per cent assurance. 
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Expected Cost of on Hand Inventory E(C) 


as a Function of the Required Number k of 
Unit Systems Up with 90 Per Cent Assurance 
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V. OBSERVATIONS AND CONCLUSIONS 


A. OBSERVATIONS 


1. Marginal Analysis may Produce only a Near-Optimal 
Se luiaon 


The marginal analysis procedure presented in 
Section III sometimes overshoots the optimum solution. 
That is, the r vector found using the procedure sometimes 
results in a higher than required availability and a higher 
expected cost of on hand inventory than needed to satisfy 
the const##aint. This behavior of marginal Binalyc#emwas 
observed during the hand calculation of the solution to 
a small problem with four unit systems and four part types 
where the constraint was E(Z) 2 3.6. Marginal analysis 
yielded an E(Z) of 3.709 and an E({(C) of $53.87. it 
was possible, however, on the last iteration to achieve 
an E(Z) of 3.601 at a cost of $48.63 by increasing a 
different a than the one associated with the maximum 
nae A dynamic programming procedure would overcome this 
shortcoming, but dynamic programming rapidly becomes compu- 
tationally cumbersome as the number of decision variables 
increases. In the sample problems discussed in Section IV 
there were 159 decision variables, and in many possible 
applications of the model the number would be in the 
thousands. Thus the simple marginal analysis procedure may 
be more computationally feasible than dynamic programming 


for many applications of interest. Further, the relative 
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differences between the optimal solution and the marginal 
analysis solution tends to decrease as the size of the 
problem increases. 

2. Program PKSMIN is Much Faster than Program EZMIN 

Sample problems 1 and 2 in Section IV resulted in 

about the same values of E(Z) andof E(C). For problem 
1 E(Z) was 47.58 and E(C) was $418.04. For problem 2 
E(Z) was 47.46 and E(C) was $387.88. On the IBM 360 
program EZMIN required over four minutes to compute the 
results for problem 1, while program PKSMIN solved problem 
2 in nineteen seconds. The reason for this difference is 
that E(Z) amore complex function than Pr(Z >k). For 
a given r vector E(Z) is the suwever k ‘from one to 
S of Pr(Z >k). A comparison of the JMAX segments of 
subroutines OPTIMR and OPTMZ shows that OPTIMR requires 
more computations at each iteration (See Flow charts B4a 


and C3a.) 


3. Functional Relationships of Expected Cost Versus 
Reqiired Availability are Readily Obtained 


The functional relationships of expected cost of 
on hand inventory versus availability discussed in Section 
IV were obtained by simply adding a DO loop to the two 
Main programs. Thus for E(C) versus A, the value of 
A. was initially set to 0.89, and at each iteration of the 
DO loop A, was increased by 0.01 for ten iterations. 
Similarly for E(C) versus k at the 90 per cent assurance 
level k was initially set at 44 and incremented by one 


at each iterations for six iterations. 


GD 


Using the functional relationships of E(C) versus. 
availability we can estimate the maximum availability which 
can be achieved if there is a constraint on E(C). Suppose, 
for example, we desired to maximize A, subjects tom EC) 

4 $2000.00. From@able yw "inySe Ctlonmieweuseeethat the 
maximum A. would be between .96 and .97. To determine 
maximum A. more precisely and to determine the associated 
r vector we could run program EZMIN for a series of values 
of Ao: says .961, —Q62y — if 39.60). euThemmes wien oebalie) 

for one of these A. - values should be very close to 
$1000.00. 

Computing time in developing the E(C) versus avail- 
ability functions can be conserved if at each iteration 
after the first the starting values of the — are the 
final values computed on the previous iteration. This can 
be ‘accomplished if. the -initializing subroutines INITAL and 
INTLZ are skipped on the second and subsequent iterations. 
This technique reduces the computing time to about half 
that needed if the initializing subroutines are called at 


each iteration. 


B. CONCLUSIONS 
1. Advantages of the Model 
The model provides the first technique of which 
the author is aware for explicitly taking into account 
system availability in the determination of continuous 
review inventory policies for repair parts. The model 


allows for variable safety levels and yields a set policies 
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where the ratio of marginal availability to marginal cost 
of inventory investment is approximately the same for all 
part types. The result is that the model tends to favor 
low-cost, high-demand items with relatively high safety 
levels. High-cost, low-demand items are less well protected. 
All items are, however, well enough protected to insure, 
under the assumptions of the model, that the availability 
requirement will be met. 

Compared to models which might more realistically 
represent the "real world" the one presented here possesses 
the advantage of computational feasibility. Programs 
EXMIN and PKSMIN produce solutions in a reasonable amount 
of time. Refinements in the programming could probably 
improve run times. 

2. Limitations of the Model 

The limitations stem primarily from the assumptions 
which were made for the sake of mathematical simplicity. 
For example, the constant lead time assumption is not 
usually true. Lead times can have considerable variability. 
Lead time variability will tend to decrease the actual 
availabilities achieved compared to those predicted by the 
model. On the other hand the model does not consider 
the fact that requisitions for parts which are causing a 
unit system to be down have a higher priority and thus a 
shorter lead time than normal replenishment requisitions. 
This fact will tend to make achieved availabilities higher 


than those predicted. 
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The optimal cannibalization assumption is seldom 
completely true. The significance of this assumption 
decreases, however, as the required availability increases, 
Since at higher availability levels we expect fewer back- 
erders and therefore fewer occasions fon, cannibalization. 

The significance of the assumption of zero or one states 
for unit systems will depend upon the stringency of the 
serviceability criteria used to determine whether a unit 
system is up or down. The more stringent the criteria, 
the more parts will have bs equal to 2° Parts with 
2s equal to 4 will generally be better protected than 
parts with bs less than aks The most stringent criteria 
possible would be that every part on the unit system must 
be working if it is to be counted as being up. Thus we 
can see that the inventory policies computed by this 
model and the resulting expected costs of on hand inventory 
are highly dependent upon the serviceability criteria for 
the unit system. 

The assumption that parts are replaced immediately with 
zero lead time makes it necessary to correct the avail- 
abilities computed by the model for down time due to finite 
time to repair. 

3. Uses of the Model 

The model was developed for use in determining 
direct support inventory policies for repaire parts for a 


given population of a particular major item. 
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Another use which suggests itself is determining 
optimal maintenance floats. Suppose it were required that 
at least k unit systems be up with probability .90. 
Assume that the capability of the maintenance unit to make 
repairs, given that needed parts are available,is fixed. 
Then the probability that k unit systems are up is 
determined by the total number S of unit systems and 
the inventory policies for repair parts. The quantity 
S - k is called the maintenance float. An interesting 
problem is determining the optimal level of the maintenance 
float. The model presented in this thesis could be used to 
estimate the optimal float level as follows. Compute the 
expected cost of on hand inventory needed to assure with 
probability .90 that k unit systems are up with zero, one, 
two, ... unit systems in the maintenance float. Then choose 
that level of maintenance float which minimizes the sum of 
the cost of the float plus the on hand inventory. 

4. Extensions 

A number of possible extensions of the model 

appear to be worthy of investigation. For example, relax- 
ation of the assumption that all unit systems are identical 
would be useful. A next higher level of complexity of the 
system would be one in which there are two types of unit 
systems. For example, one type of unit system might be a 
particular type of weapon. The other unit system might be 
a piece of fire control equipment which controls the fires 


of several weapons. If a fire control unit is down all the 
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weapons it controls are also down. Application of the 
approach of this node to optimal inventory policies for 
such a system seems to offer interesting possibilities. 

In any case it is hoped that the model presented here 
will prove useful in the development of system availability- 
oriented operating policies for military repair parts 


inventories. 
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APPEND LE a 
CONVENXITY CONSIDERATIONS 


In ogder to insure the Ooptimaletysct the & _ voever 
found using the procedure outlined in Section III it is 
necessary that the cost function be convex and that the 
region in r-space defined by the constraint be convex. 

rine Cost Function 

The objective function in the model presented in 
this paper is the expected cost of on hand inventory, 


which is given by 


D(C) = 


Ine 


CD (©. 
gP (Q50F5) 


j=l 


Since the Q. are fixed by the operating level of supply 
policy as described in Section III, we are concerned 

only with convexity with respect to the eae f=" hess. Fee 
Further, since E(C) is a sum of terms each of which is 

a function only of ras E(C) is convex if each D,(Q. 5) 
1s convex. Since Q. is fixed let us drop it as an argu- 


ment; also let us fix Jj and drop the subscript. Hadley 


and Whitin [Ref. 8, p. 184] show that 
Dia) = (Qete l)/2 + «x + Bye) 


where B(Q,r) is the expected backorders, given a_ (Q,r) 


policy. D(r) is linear in r except for the B(Q,r) term. 


a 


Thus D(r) is convex if B(Qye)e 1SMeenveswin® r.> Agcin, 
since Q is fixed, let us drop it as an argument. Hadley 
and Whitins [Ref.n Sy. p. b8at shew Vea 


co 


B(r) = (1/Q) % ylP(ytr+l;u) - Ply+r+O+1;1)). 


y=0 
cee 
AB (r) = B(r+l) - B(x), 
and 
A“B (r) =e Beer tie) See Emteie) 
i i i 


Dit) AtSeconvexa ur A<B (r) < 0 for all xr. Since the par- 


—= 


ameter wu is fixed for fixed Jj, let us further simplify 


the notation by letting 


P(x) = Peron; 


pS) eS ep aa) 


Further, let 
Ee (22) = ©8 (is) . 


Mien” D(r)iPasemeonvex for ree r* if A°B* (x) ~ @ + formall 


1g es et Nel) 


= B(ytmtZ) SP (ye ez) 
0 -P(y+rt+1) + P(yt+r+Qt+1) 


A B* (x) 


8 


= Y yliplytr+OQt1l) - p(ytrt+l)]. 
Vou 
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Note that we Have dropped the term Mee y = 0 since Ft 


contributes nothing to the sum. 


Now 
2ae 7 - pPivytr+O@e2) —sp (ye) 
aes Hii. Y)-ply+r+O+1) + P(ytr+1) 
= Z£ yflAply+r+O+1) ="“Ap(ytrmrly] 
ye 
where 
Noi = pel) =p (x) 
F y OFT) a , i oe 
- (xcttaly): ! x! 
eV aaa 
~ =e RIES 
Thus 


Z _ ia Ls VseOm2 jie Gai2 
ASB (as) sae _ y [ Tron p (y+r+Qt1) ytrto p(y+rt+1) ] 


y=l Vata 7 iar 2 
(r+3-U) _ $B Ra.) 
a ia ls “=rege3 PERF?) 
+ bere (re3) - 2SUERORENY 5 (r+g43) 
+ (Q+1) AEFRF SI» (r+9+2) Te ee 
“= 
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Os ce 
=, y (TEXTE) p(ytr+1) an ae (FRYE SEE) p (rty+1). 
y=1 a fal Rd 
Now every term in this series is non-negative if r > u-3. 
Thus D(r), is convex in r for Yr =e ee buethesanore. 
numerical calculations show that for values of u up to 
40 that if Qo> u-.5 then D(r) is convex for r> -l. 
This numerical result is probably due to the fact that the 
mode of the Poisson probability mass function is approxi- 
mately u-l, which means that the most negative values 
of (rt+yt2-u)/(rt+y+2) in the series are multiplied by 
small values of p(x;u). When r 2-1, all terms in the 
series for y >u- 1 will be positive. Thus it is 
likely that the sum of the positive terms in the series 
is greater than the absolute value of the sum of the cout 
for which y +1- u is negative. 
Thus, as far as the cost function is concerned, we can 
be sure that the @tunetion Is convex in) «forearm! 
YE Cee «=Wiere™™ a= Uprese rye and we have reason to 
believe the function is convex for r2> -—-I1 when OQ; 2 Use 
j= lope 
2. The Constraint Function 

The two alternate forms of the constraint functions 
are, first acrequired expected number of operational unit 
systems, expressed by 

S N (Sa;-kb.) 


5, ala ee ae 
k=1 j=1 y,=0 jv4 . 
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and second a required probability that at least k unit 


systems are operational, expressed by 


N Saj-kb. 
II »: Yo ly> ee) P nin ‘*) 
Jao a 0 enn 


If the constraint region is to be convex, the above func- 


tions must be concave functions of r= (Ly pee rd Note 


N° 
that the left hand side of the first constraint is the sum 
over k of terms of the same form as the left hand side 
of the second constraint. Thus if we can show the con- 
ditions for concavity of the left hand side of the second 
constraint for a general k, we can easily show the 
conditions for concavity of the first constraint. 

The difficulty in showing the concavity of the second 
form of constraint is that it is not a separable function 
of the Se It can be transformed into a separable function 


by taking the logarithm, in which case it becomes a sum of 


terms, each of which is a function of only one Lae Iwie.; 


N 
InfPr{z# > kK}) = * injeee 2 
ia : 
where 
Sa4-kb. 
Pr{Z. > k} = 12 2 (ya) 
a | es 
so 


ne 


Let 


" N 
{eye = Pr{Z > heal 


fo aea)e, 
: Ae 
; y 


1 


£.(r.) = Pri4Z] eek tee 
7 7. { 52 re 


é 
C] 


d (r) = In ioe, 
a sb a = ln fe) 


Now it can easily be shown that if eco 1S concave, 
aL is also concave, and thus $(r) is concave. Un- 
PoreEunace ly, the Concayily. Clad) a) OCS slo tad msn aan 


concavityget f(r) .. Noweconordcr 


Sa ~kbs 
Fala) a ee) Y ioe Vantaa sania) 
y.=0 ) aa 
J 
Let 
M= Sa. - kb., 
o 
then 


EgeGee) = Yai > 7.) eee te OM ee) 
a ;) 7 7 5 { =) 


In section III we showed that 


Thus we have 


M 
fGen) = ¥. s-l1) + 
5% %? E 5 ae a 


Now let 


Maes. ie.) = £2, ireiee— fe at 
pt, le) = 3 lee oie 


J 
then 
Mtr F2 Mt+r.+1 
a ee: w. (ea =o pw Aaa tae a} 
ry J u=0 J u=0 J 
= Y,(M+tr.+2;-1). 
J J 
Let 
Re F. = ANEW) — A 
1 en J Ge J 
J J J 
then 
Ne fF. = YY. (Mtr+3;-1) - VY. (Mtrt+23;-1) 
rs. J J J 
J 
P(Mt+r4+3) = P(M+Q5+r.+3) 
= (1/0) J 
-P(M+r+2) - P(M+Q+r+2) 


(1/Q) [p (M+Qt+r+2) - p(Mtr+2) J 


Now t, is concave for all re a ae if the expression 


in brackets above is < 0 for all 4, ES eae Thus we seek 


FS suieh that for all t. 28s 


a 


Mtr.+2) > M+O.+r.+2 
p ( ee ) 2. pK 2. ee ) 


or 
Mt+r.+2 M+Q.+r.+2 
‘igen 1 ) vu i: Q. be ) cu 
J... s 2 
Mtr .+2)! a M+Q.+xr.4+2)! 
( : ) Q. 5 ) 
or 


(Mtr .+Q.+2)! 
< 
a BESS ach) 


The right hand side of the above inequality is an increasing 


Pe © hel) esses O., and M. It is usually true that the 


1 =? 
operating level policy used within the Army is at least as 


many days of supply as the procurement lead time. This 


means that Q. aa ee Assume this is true. Now if k= S§ 


and Be = ae then M= 0. If eer 1S concave, 
a! will also be concave, Since if the inequality 
holds for M= 0, it will hold for M>0. Let M = 0. “Gites 


the inequality becomes 


(r.+O7g2oe! 
ae cha 
a (r.+2)! 
J 
For example, if ee = - 1 and Q. = 1, then ts must be 
equal to or less than 2. In other words, for us <oa2 | se 
1S concave in ie for all a >- il, Q. > 1, and k = 3s8 


To insure the concavity of Sd For iu.> 2, in thie 


FORTRAN program for the E(Z) = A.S form of the conetramme, 


78 


the subroutine INITAL finds an initial set of ty such that 
the concavity conditions for k= S are satisfied for — 

equal to or greater than the initial og Similarly, in 
the FORTRAN program for the Pr(Z=k) = ea form of the 


constraint, subroutine INTLZ finds an arias 1 ae) set of te 
such that each pr(Z. = k; r,) is a concave function of 
= Fox a equal to or greater than the initial value. 

The solution procedure employed in the computer programs 
thus insures that each £5 (rs) 1S concave in x in the 
region that the marginal analysis procedure searches for 
an optimum. The question of the concavity of E(e)= Ml f(xy 
remains open. To date the numerical results have ae inid- 
cated that there 1S a problem, i.e., the numerical results 


are quite reasonable. No proof, however, of the concavity 


OL rar) has been found. 


oS 


APPENDIX 5 
FLOW CHARTS FOR FORTRAN PROGRAM EZMIN 


FORTRAN program EZMIN was developed to solve the problem 
of finding the vector r which minimizes the expected cost 
of on hand inventory subject to a required expected number 
of unit systems up. This appendix contains flow charts for 
the main program and the following subprograms: 

Subroutine PSITAB, which computes a table of values of 


Pm ,Y) forme je= Gee, .,Neeeand Y = 1,2, LOOM where 


Pia , Ys ¥5 (Y-lixr, = =— 9) 


Subroutine INITAL, which computes a set of initial 
values of ae such that cae ~ ae. ge) Slee ee 
Subroutine OPTIMR, which uses a marginal analysis tech- 


neque to find an oOpemmal “vector, 
Function PRZJK, which computes Priz. = ks ryt. 
Function EBO, which computes af v4). 
Function PPT, which computes individual terms of the 


Poisson probability distribution. 


Function PCUMT, which computes complementary cumulative 


terms of the Poisson probability distribution. 
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NO .@OF UNIT SYSTEMS 
NO..#: GF PART nis 
AZ = REQUIRED AVAILABILITY 
OL = OPERATING LEVEL 
ALEVEL = ACTIVIT, GEVED 


COMPUTES TABLE OF PI(j,k) = j(k-lir, = - 1) 
FOR Jj = ~am@e Ne k = Ly... , 200 





Figure Bl. Flow Chart for EZMIN Main Program 
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pba, = AZ*S AND Pr (2; = Sirs) IS CONCAVE IN LS 






GPT IMR 


COTE Oana le emma oy eer serene Ne) See le etree 
EAC) IS MINAMIZER BUBIJECTSTO. E(Z) = Azee 


a OP 


82 


START 









KMAX(J) = 100 


eo. a 30 - 
| a 1,100 


K >KMAX (J>——— 


= Q. 


| 

| 

| 

| 

| 

| Padi) 9 =) OMe ie Oe a) 
| nr 
| 

| 

| 

| 

| 










K >u, 
amcde Pua( J ,K) 


es 


RETURN 


Figure B2. Flow Chart for Subroutine PSITAB 
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PL(J oy = 1. -. SUMO. 


yes 





PI(J,K) = lie coe = B (Rone 05 t.) 1/08 






es 


KMAX (J) = K 


Figure B2. Continued 
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START 









P, = Po 


Pr (Z.=S;r.+3) 


J 
ae et 
J J 


Figure B3. Flow Chart for Subroutine INITAL 
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N 
COMPUTE Pr(Z>k) II 


Pr Oe ir.) 
1 J J 
FOR k 





Eg Z ) 


iI 
lun 
ry 
Ky 
N 
IV 
~ 


COMPUTE a = EE ae 


il 
ios 8 


- 


COMPUTE prCge™. = eee to BROKE) = ae J = 1,.. Nn 


ITER = 0 


yes € 
7a a) 


Pak = T°. 


Find jmax = {j*: we = max Aas 


E(Z;r5+1) ~ a a) 
where A. = DTCOST 
and c Be (2. >k3r2+1) 
Bis oad a pot eee er (Zn 2k; Xr.) 
J J 
Seb. £(4) = E(a a + 1) 
jmax 


(See Figure B4a for details of this block) 


Figure B4. Flow Chart for Subroutine OPTIMR 
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Jes Ts + J 
jmax jmax 
ED max = EBOr 


ma ee 


jmax) 


(1 + EBO(r 


jmax jmax 


COMPTE EOH., = +) oe ee - Ue Ney = 
5 [(Q. )/ 3 5 4! 3 


WROTE: Ie (2), EB(C), 





RETURN 


Figure B4. Continued 
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oe 


3 
ae —_—_ —_—=_ — —— ——— an —_— tS atestiaaaell ——= 


PKNEW = Pr(z>K) 





— ee eee eee 


DELTAJ = (EZTRY - ta@Z) ) /DT COST: 
yes ma > 
no 
rf 
>DHOLD 
ves 


DELTAJ 


EAiRY 





E(Z) = EZSAVE 


Figure B4a. Flow Chart for JMAX Segment of OPTIMR 


88 






START 











KHa f° OA. - kb. + ri 
J J J 


yes 








PRZJK = PRZJK + °PT (J) RY 


+ a © 
RETURN 


Figure B5. Flow Chart for Function PRZI5K 
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EBO = 0. 
= RETURN 
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Y= 0 





Y+e4 
K- = Ser ts + 2 
EBO = EBO + Y*PI(J,K) 








yes 


Figure B6. Flow Chart for Function EBO 


START 


MU = RHO*T 


eS 


PPT = DEXP(-MU) 
no 


RETURN PPT = DEXP (-MU) 


<a> 
| 


PPT = PPT*MU/DFLOAT (TI) 





La a84@ 
RETU RN 


Figure B7. Flow Chart for Function PPT 


el 


START 


MU = £RHO*T 


PCUMT = 1.D00 PCUMT = 1.D00 - 
DEXP (-MU) 
RETURN 





SUM = DEXP (-MU) 
PROD =sc50R 
THI = IX-1l 


PROD = PROD*MU/DFLOAT (T) 
SUM = SUM + PROD 





PCUMET = 1.800, 2550UM 
RETURN 


Figure B8. Flow Chart for Function PCUMT 


APPENDIX C. 
FLOW CHARTS FOR FORTRAN PROGRAM PKSMIN 


FORTRAN program PKSMIN was developed to solve the prob- 
lem of finding the r vector which minimizes the expected 
cost of on hand inventory subject to a required probability 
that at least k unit systems are up. This appendix con- 
tains flow charts for the main program and the following 
subroutines: 

Subroutine INTLZ, which compute a set of initial values 


of r5 such that as ak) > PG > K) j=l,.. 


Wy? 
Subroutine OPTMZ, which uses a marginal analysis tech- 
nigue to find an optimal _.r vecter. 
Program PKSMIN also uses subroutine PSITAB and function 
subprograms PRZJK, EBO, PPT, and PCUMT, flow charts for 


which are in Appendix B. 
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START 


READ: .S. = NO. OF Uhm SYSmErvs 
KS = REQUIRED NO. OPERATIONAL 
UNTT SYSTEMS 
N = NO. OF PART TYPES 
PKSMIN = REQUIRED ASSURANCE 
OL = OPERATING LEVEL 
ALEVEL = ACTIVITY LEVEL 


WRITE: 5S, KS, N, PKSMIN, OL, ALEVE 


COMPUTES TABLE OF PI (J,k) 
FOR 47s...) ieee 





Figure Cl. Flow Chart for PKSMIN Main Program 
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OPTMZ 
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Figure Cl. Continued 
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Figure C2. Flow Chart for Subroutine INTLZ 
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START 


N 
COMPUTE P(2>KS) = I 
j=l 











Pr (Z.>KS) 






GOMPUTE EB. = EBO{r.) =. yee) 5 — ue 
j j jee 3 j 
y 
COMPUTE DTCOST; = c,(1 + EBO(r +1) - EB,),j = 1,...,N 
ITER = 0 
“ ey 
i (8) 


ITER = LEER?) 


jmax = {j*: 4." = Be Aa). 
Pr (Z>KS;r.+h) - P(Z>KS) 
Cine : 
Pr(2i2KS;r.+1) 


J 
Prw.or . i. 
ae) 


rae ree =—P ria ee) 


set Pr (Z>KS) = Pri ZoKos oon 


jmax sa) 


(See Figure C3a for details of this block) 





Figure C3. Flow Chart for Subroutine OPTMZ 
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jmax jmax . 
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RETURN 


Figure C3. Continued 
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M= Sa. -— KS*b. + t = |S 


es PKSTRY = Pr(Z>KS) 


PKSTRY = Pr (Z>KS) [Pr (2,>KS)+P1(J,M) ]/Pr(Z,>KS) 
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os Get 
no 


DELTAJT no 
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ew 
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